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Abstract — Future communication networks are expected to 
feature autonomic (or self-organizing) mechanisms to ease de- 
ployment (self-configuration), tune parameters automatically 
(self-optimization) and repair the network (self-healing). Self- 
organizing mechanisms have been designed as stand-alone en- 
tities, even though multiple mechanisms will run in parallel 
in operational networks. An efficient coordination mechanism 
will be the major enabler for large scale deployment of self- 
organizing networks. We model self-organizing mechanisms as 
control loops, and study the conditions for stability when running 
control loops in parallel. Based on control theory and Lyapunov 
stability, we propose a coordination mechanism to stabilize the 
system, which can be implemented in a distributed fashion. The 
mechanism remains valid in the presence of measurement noise 
via stochastic approximation. Instability and coordination in the 
context of wireless networks are illustrated with two examples 
and the influence of network geometry is investigated. We are 
essentially concerned with linear systems, and the applicability 
of our results for non-linear systems is discussed. 

Index Terms — Self-Organizing Networks;Self- 

Organization;Self-Optimization;Self-configuration;Stochastic 
approximation ; Stability Coordination ; Dynamical Systems ; 

I. Introduction 

Deployment, optimization and maintenance of communica- 
tion networks are complex tasks which occupy thousands of 
network engineers everyday. In order to ease this burden and 
reduce costs of operations, researchers and industrials have 
proposed to include autonomic functionalities in future net- 
works. Networks with autonomic entities are often called Self- 
organizing networks (SON). Self-organization should enable at 
least partial automation of the configuration of newly deployed 
network node (self-configuration), of parameter tuning (self- 
optimization) and of reparation of faulty network nodes (self- 
healing). 

In wireless networks, SON functionalities such as Inter-Cell 
Interference Coordination (ICIC), load balancing, management 
of Base Station (BS) sleep mode, mobility management and 
drive test automation have been identified as important use 
cases (HI)- Future standards for wireless networks such as 
Long Term Evolution (LTE) feature SON functionalities. In 
previous contributions, SON functionalities have been de- 
signed as stand-alone entities (0, 0, 0, 0, 0). Such 
contributions did not take into account the interaction between 



different SON functionalities operating simultaneously. Since 
SON functionalities are numerous, a large number of SON 
entities might operate simultaneously in the network. This 
poses a fundamental stability issue, if we cannot guarantee 
stable interaction of multiple SON entities. It is fair to say that 
a robust coordination mechanism will be the main enabler for 
large-scale deployment of SON. 

Our contribution: The main contribution of this paper is to 
propose a generic mathematical model for the interaction of 
multiple SON mechanisms running in parallel, along with a 
practically implementable coordination mechanism. We model 
SON mechanisms as control loops, so that the system can 
be described by an Ordinary Differential Equation (ODE). 
Stability is studied using the Lyapunov approach used by con- 
trol theorists. In particular, stability conditions can be stated 
in terms of linear matrix inequalities (Q). When stability 
does not hold naturally, we propose a coordination mechanism 
which forces stability. The design of a coordination mechanism 
corresponds, in the context of control theory, to the concepts 
of controllability and state-feedback synthesis. We propose a 
generic coordination mechanism which can be implemented 
in a distributed way. When the network is controlled using 
measurements corrupted by additive noise, the results remain 
valid using stochastic approximation theorems (0). 

To the best of our knowledge, no previous contribution 
had studied the problem of SON coordination using this 
control theory/stochastic approximation-based framework and 
provided a generic coordination mechanism which is practi- 
cally implementable (distributed and possibly asynchronous). 

The rest of the paper is organized as follows: in Section Ull 
we state the proposed model for interaction of SON mech- 
anisms running in parallel and the coordination problem to 
be solved. In Section [HI] we examine the case where perfor- 
mance indicators are linear functions of the parameters, and 
propose a practically implementable coordination mechanism. 
In Section [IV] we illustrate the application of our model to 
traffic management in wireless networks with two examples. 
Section [V] concludes the paper. In appendices [A] and [B] we 
recall the basic notions of stability for ODEs and linear ODEs 
respectively. 
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II. SON COORDINATION AS PARALLEL CONTROL LOOPS 

A. The model 

A SON mechanism is an entity which monitors a given 
performance indicator and controls a scalar parameter. The 
current value of the performance indicator is observed, and 
the parameter is modified accordingly to attain some ob- 
jective. We consider I > 1 SON mechanisms operating 
simultaneously. We define 0, the parameter controlled by the 
i-th SON mechanism and 9 = (0i,...,0j) the vector of 
parameters. The i-th SON mechanism monitors a performance 
indicator and updates its parameter 0, proportionally to Fi{9). 
F{9) = (Fi(0), . . . , F 7 (0)) is the direction of update of 9. 

We say that the i-th SON mechanism operates in stand- 
alone mode if all parameters but 9i are kept constant. Namely 
all the other mechanisms are shut down. The i-th SON 
mechanism operating in stand-alone is described by the ODE: 

e\ = Fi{0), 

o../ •• <•■ (i) 

We say that the SON mechanisms operate in parallel mode if 
all parameters are modified simultaneously, which is described 
by the ODE: 

9 = F(9). (2) 

We say that the i-th SON mechanism is stable in stand-alone 
mode if there exists 9*' 1 which is asymptotically stable for (Q]i. 
The definition of asymptotic stability is recalled in appendixlAl 
It is noted that 9*' 1 can depend on 9j, j ^ i. We say that the 
SON mechanisms are stable in parallel mode if there exists 
9* which is asymptotically stable for (01. Typically, the SON 
mechanisms are designed and studied in a stand-alone manner, 
so that stand-alone stability is verified. 

However, stand-alone stability does not imply parallel sta- 
bility. First consider a case where Fi(9) does not depend 
on 9j, for all pairs (i,j) , j ^ i. Then (ffji is a set of 
/ parallel independent ODEs, so that stand-alone stability 
implies parallel stability. On the other hand, if there exists 
i 7^ j such that Fi(9) depends on 9j, then the situation is 
not so clear-cut. We say that SON i and j interact. Namely, 
interaction potentially introduces instablity. 

In the remainder of this article we will be concerned with 
conditions for parallel stability, and designing coordination 
mechanisms to force stability whenever possible. 

B. Examples of parallel mechanisms 

Two particular cases of parallel mechanisms will be of 
interest. The first case is what we will call zero-finding 
algorithms. Each SON mechanism monitors the value of a 
performance indicator and tries to achieve a fixed target value 
for this performance indicator. Namely: 

Fi{6) = M9) -/i, O) 

where /j is the performance indicator monitored by SON i and 
/j - a target level for this performance indicator. The goal of 
the t-th SON mechanism is to find 9* such that fi{9*) = f i . 



If 9i i-> fi(@i, . . . ,9i, . . . , 9i) is strictly decreasing 1 < i < I 
then stand-alone stability is assured. 

Another case of interest is maximization algorithms. Each 
SON mechanism tries to maximize a given performance in- 
dicator. There exists a continuously differentiable function gi 
such that: 

F l {9) = V ei5i (0). (4) 

In stand-alone mode, SON i indeed converges to a local 
maximum of gi. If we restrict 9 to a closed, convex and 
bounded set and if 9i n- gi{6\, . . . , 6%, . . . , 0j) is concave 
1 < i < I , we fall within the framework of concave games 
considered in |9). An important result of (9| is that if we add 
an assumption called diagonal strict convexity, then parallel 
stability occurs. However, diagonal strict convexity is fairly 
restrictive, and without it there is no guarantee that parallel 
stability occurs, and coordination is needed. 

C. Discrete time algorithms 

Our model is based on ODEs, which are both determin- 
istic and continuous-time objects. In a practical system, the 
parameters evolve in discrete time. Furthermore, the param- 
eters evolve stochastically because they are updated based 
on noisy feedback, due to measurement noise. Denote by 
9[t] = (0i [t], . . . , 9i[t]) the value of parameters at time 
t G N. At each parameter update, performance indicators are 
estimated based on measurements so that SON i obtains the 
quantity Fi(9[t\) + Mi[t], with E[Af,[t]] = 0. The additive 
noise {Af[i]} tgN is introduced by the measurements. The 
parameters are updated using the noisy feedback: 

9[t + l]=6[t]+e(F(6[t])+M[t)). (5) 

where e > is a small constant. Stochastic approximation ([8], 
1 10]) gives a strong link between the noisy, discrete-time © 
and the ODE studied in our model ©. In particular it can 
be shown that (under some technical conditions) {0[£]} tgN 
converges to asymptotically stable sets of the ODE when 
e — > + . The main point is that our model based on ODEs 
is sufficient to study the stability of the discrete-time, noisy 
algorithms used in practical systems. 

III. Coordination for linear systems 

A. Linear systems 

In most of the remainder of this paper, we will study the 
case where F is affine: 

F{9) =A9 + b, (6) 

with b a vector of size / and A a matrix of size I x I. We 
assume that A is invertible and we define 9* = — A~ 1 b. The 
SON mechanisms running in parallel are described by the 
linear ODE: 

9 = A9 + b = A{9-9*). (7) 

It is noted that in the linear case, we always fall within 
the scope of zero-finding algorithms described previously, by 
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defining: 



M9)= A ^ 
i<j<i 

7 1 = -k, 

0i = M6)-7 i . 



(8) 

(9) 
(10) 



In particular, stand-alone stability occurs iif < , 1 < i < 
I, i.e all the diagonal terms of A are strictly negative. Basic 
results on linear ODEs are recalled in appendix [B] Namely, 
parallel stability holds iif all the eigenvalues of A have strictly 
negative real part. 

B. Coordination 

1) Coordination mechanism: If A has at least one eigen- 
value with positive or null real part, convergence to 9* does not 
occur, and a coordination mechanism is needed. We consider 
a linear coordination mechanism, where A is replaced by CA 
with C a / x I real matrix. The ODE for the coordinated 
system is: 

6 = CA(6 -6*). (11) 



Define a as: 



i<3<! 



(12) 



The coordination mechanism can be interpreted as transform- 
ing the performance indicator monitored by SON i from /, to 
a linear combination of the performance indicators monitored 
by all the SON mechanisms Cj. As explained in appendix IB~1 
stability is achieved if there exists a symmetric matrix X such 
that: 



(CA) T X + XCA ~< 0, 
-< X, 

where -< denotes positive definitiveness for symmetric 
ces. In particular, 

V(9) = (e-6*) T X(6-6*), 



(13) 
(14) 

matri- 



(15) 



acts as a Lyapunov function. 

2) Distributed implementation: It is noted that the choice 
for the coordination matrix C is not unique. For instance 
C = A^ 1 ensures stability. For the coordination mechanism to 
be scalable with respect to the number of SONs, C should be 
chosen to allow distributed implementation. We say that SON 
j is a neighbor of SON i if ^ 0. We define Zj the set 
of neighbors of i. The coordination mechanism is distributed 
if each SON needs only to exchange information with its 
neighbors. 

We give an example of a coordination mechanism which can 
always be distributed. The mechanism is based on a separable 
Lyapunov function. Define the weighted square error: 



V(6) = J2Mfi(0)-fi 



r) T A T WA{6 - 9*), 
(16) 



with {wi} 1<i<:l strictly positive weights and W the diagonal 
matrix witrT diagonal elements {wi} 1<i<:i . Coordination is 
achieved by following the gradient of —V so that V is a 
Lyapunov function: 

9 = -V e V(6) = -A T WA(6-6*). (17) 

Namely, we choose C = —A T W. We can verify that the 
mechanism is distributed: 



d6i 



/i) = £a^C/iW-/i). 

(18) 

Indeed, the update of 0i only requires knowledge of ^f- and 
fj(0) — fj, for j and this information is available from 
the neighbors of i. 

It is also noted that such a mechanism can be implemented 
in an asynchronous manner, i.e the components of 6 are up- 
dated in a round-robin fashion, or at random instants, and the 
average frequency of update is the same for all components. 
The reader can refer to ifTTI [chapters 6-8] for the round- 
robin updates and [8][chapter 12] for the random updates. 
Asynchronous implementation is important in practice because 
if the SONs are not co-located, maintaining clock synchroniza- 
tion among the SONs would generate a considerable amount 
of overhead. 

C. Applicability of the linear model and parameter estimation 

For practical systems, performance indicators F{&) are not 
linear functions of 6. However, as long as they are smooth, 
they can be approximated by linear functions using a Taylor 
expansion. Consider a point 6* with F(6*) = 0. If the values 
of are restricted to a small neigborhood of 6*: 



F{9) w JF(0*)(0 - 9*), 



(19) 



with JF(9*) the Jacobian of F evaluated at 9*. The Hartman- 
Grossman theorem (pl2|) states that on a neigbourhood of 9*, 
stability of the system with linear approximation implies sta- 
bility of the original, non-linear system. Hence implementing 
the proposed coordination mechanism where A is replaced by 
JF ensures stability if we constrain 9 to a small neighborhood 
of 9*. 

The parameters A and b might be unknown, and we can 
only access to noisy values of F(9) for different values of 9. 
The crudest approach is to estimate A and b through finite 
differences: 

_ fiifi + eiSBi) - fj{fi - eiSBi) 



289, 



h w /i(o). 



(20) 
(21) 



with e; the z-th unit vector. The results are averaged over 
several successive measurements and additive noise is omitted 
for notation clarity. In general, the measurements of F are 
obtained by calculating the time average of some output of 
the network during a relatively long time, so that a form 
of the central limit theorem applies and the additive noise 
is Gaussian. In this case, a better method is to employ 
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least-squares regression. Least-squares regression is a well 
studied topic with very efficient numerical methods ( IT31 ) 
even for large data sets so that estimation of A and b is not 
computationally difficult. 

Finally, since practical systems do not remain stationary for 
an infinite amount of time, a database with values of A and 
b for each set of operating conditions must be maintained. 
In the context of wireless networks, the relationship between 
parameters and performance indicators changes when the 
traffic intensity changes because of daily traffic patterns. For 
instance, during the night traffic is very low, and traffic peaks 
are observed at the end of the day. A database with estimated 
values of A and b at (for instance) each hour of the day should 
be constructed. 

IV. Application to wireless networks 

In this section we illustrate instability and coordination in 
the context of wireless networks using two examples. We show 
that instability occurs even in simple models with as few as 
two SONs in parallel. 

A. Admission control and resource allocation 

1) Model: We consider a BS in downlink, serving elastic 
traffic. Users enter the network according to a Poisson process 
with arrival rate A, to download a file of exponential size 
a, with E [a] < +oc. The BS has x max parallel resources 
available, and we write x G [0, x max \ the amount of resources 
used. We ignore the granularity of resources, either assuming 
that there are a large number of resources or using time 
sharing, using each resource a proportion x of the time. 
Depending on the access technology, resources can be: codes 
in Code Division Multiple Access (CDMA), time slots in Time 
Division Multiple Access (TDMA), time-frequency blocks in 
Orthogonal Frequency-Division Multiple Access (OFDMA) 
etc. When a user is alone in the system, his data rate is Rx. 
Users are served in a processor sharing manner (for instance 
through Round Robin scheduling): if there are n active users, 
each user has a throughput of Admission control applies. 
We define b > a blocking threshold and the probability of 
accepting a new user when n users are already present in the 
system is cf>(n — b) with (j> : M — > [0, 1] a smooth, strictly 
decreasing function and <fi(n) —> 0. We choose as a 
logistic function for numerical calculations: 



conditions: 



1 



1 



(22) 



Define n(t) the number of active users in the system at time 
t, then n(t) is a continuous time Markov chain. n(t) is ergodic 
because the probability of accepting a new user goes to as 
n — > 00. We define the load: 



p(x) 



AE [a] 
xR 



(23) 



We write 7r the stationary distribution of n(t). n(t) is re- 
versible, and 7r can be derived from the detailed balance 



7r(n, x, b) 



z n > Q p(x) n n n k =o^k-by 



(24) 



Using Little's law, the mean file transfer time is given by the 
average number of active users divided by the arrival rate: 

1 x - 



T { x , b ) = J 2^ n7T ( n > x ' & )- 

n>0 



(25) 



Let R m in a minimal data rate required to ensure good Quality 
of Service (QoS). We say that there is an outage in a state of 
the system if users have a throughput lower than R m in- When 
there are n active users in the system, outage occurs if and 
only if: 

n>^. (26) 



Rr, 



The outage probability is then: 

0(x,b) = y"V(n, x, &)l(o,. 



xR 

Rmin 



(27) 



In this model, x — > 0(x, b) is not smooth, which is why we 
introduce the smoothed outage O: 



0{x,b) 



n>0 



7r(n, x, b)ip I n 



xR 

Rm.in 



(28) 



with tp a smooth function approximating l(o,+oo)- We also 
choose ip as a logistic function for numerical calculations. 

This queuing system is controlled by two mechanisms, and 
that control occurs on a time scale much slower than the 
arrivals and departures of the users, so that the mean file 
transfer time and outage probability are relevant performance 
metrics, and can be estimated from (noisy) measurements. The 
mechanisms are: 

> Resource allocation: a mechanism adjusts the amount 
of used resources to reach a target outage rate. Such 
mechanisms have been considered in green networking 
when a BS can switch off part of its resources in order 
to save energy. 

Another application is interference coordination: using 
more resources will create inter-cell interference in neigh- 
boring BSs and degrade their QoS. Hence BSs should use 
as little resources as possible, as long as their target QoS 
is met. 

> Admission control: another mechanism adjusts the ad- 
mission control threshold to reach a target file transfer 
time. In particular, it is noted that without admission 
control, the mean file transfer time becomes infinite in 
overload. 

It is noted that x — > 0(x, b) is strictly decreasing and 
b — > T(x,b) is strictly increasing. Using the notations of 
Section HU we have 1 = 2 control loops, 6\ = x, 62 = b, 
fl = O, ji = — T. Consider 9* = (x*,b*) an operating point. 
The stability in the neighborhood of (x* , b*) can be calculated. 
The system will fail to converge to the desired operating point 
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as long as the Jacobian matrix has a negative determinant, so 
that there are two eigenvalues of opposite sign: 



dOOT, . ... dOdT. . n 



(29) 



2) Results: We now evaluate the stability region of the 
system numerically by checking condition d29b for various 
operating points. We choose the following parameter values: 
A = 0.5users/s, E[a] = WMbits, R = 15Mbits/s, 
Rmin = 2Mbits/s, x ma x = 1- Figure Q] presents the results. 
In the white region the system is stable, and in the gray region 
it is unstable. Even in such a simple setting with 1 BS and 2 
SON mechanisms, instability occurs for a large set of operating 
points. 




Ressouce utilisation (%) 

Fig. 1 . Stability region of the system 

B. Distributed interference control 

1 ) Model: We turn to a case in which the SONs are not co- 
located, and distributed coordination is required. We study a 
model in which each BS can choose its transmitted power 
in order to control the QoS in the network. We consider 
a dense network, and the objective of each BS is to make 
sure that their neighboring BSs achieve a target coverage 
probability, by not transmitting at too high power. We define 
the coverage probability as the proportion of users whom 
achieve a given minimum data rate. This is relevant to current 
wireless networks, where BSs are linked with an interface 
(called X2 interface in the LTE standard), so that a BS 
experiencing low data rates or congestion can send an alarm its 
neighbors to force them to transmit at lower power to reduce 
inter-cell interference, or to offload some of its traffic to reduce 
congestion. 

We consider N s BSs serving a bounded area A. The zone 
served by BS i is written Aj. Interference from neighboring 
cells is treated as Gaussian noise. We denote by hi(r) the 
signal attenuation between BS i and location r. hi(r) includes 



both path-loss and shadowing. For numerical calculations, we 
use the classical model where the signal attenuation (in dB) 
is a linear function of the distance, and shadowing is (in dB) 
a centered Gaussian random variable. Fast-fading is ignored. 
BS i transmits at power Pi. The Signal to Interference plus 
Noise Ratio (SINR) at location r G A s is calculated by: 

hi{r)Pi 



S l {r) 



(30) 



with Nq the thermal noise. We denote by Ri(r) the data rate 
at location r when no other users are being served by BS i. 
The data rate iij(r) is calculated using the Shannon formula: 



R i (r)=wlog 2 (l + S i (r)), 



(31) 



with w the system bandwidth. Consider a target data rate 
Rmin- We define K{ the coverage probability for BS i, which 
is the probability that a user arrives in Aj at a location r 
such that Ri(r) > R m in- The coverage probability can be 
calculated by: 



l {R,(r)>R, 



iri y(r)dr. 



(32) 



Consider users arriving in the network according to a space- 
time Poisson process of intensity A, i.e the number of users 
arriving during time dt in a region of size dr centered at 
location r is \drdt. Then the number of users arriving during 
time interval [0,T] in Aj which have a data rate superior 
to Rmin, divided by |Aj|AT is an unbiased estimate of Ki. 
We define £>j the set of neighbors of BS i. The neigbouring 
relation can either come from geographical proximity, or from 
propagation conditions, so that neighboring BSs are BSs which 
strongly interfere each other. 

We define Gj the coverage probability of the neighbours of 
BS i by: 



Gi = 



(33) 



We study the case where each BS adjusts its transmitted power 
Pi to avoid degrading the network QoS and make sure that 
its neighbors reach a target coverage probability Gj. BS i 
transmitting at higher power decreases the data rate of users in 
Aj , j £ Bi, so that Pi n- Gj is strictly decreasing. It is noted 
that we work with Pj in dB (not in linear scale). Using the 
notations of Section HH we have I = N s SONs, with 9 t = Pi, 
and fi = Gi. 

2) Results: For numerical calculations, the signal attenua- 
tion (in dB) at distance d (in km) is 128 + 36.4 log 10 (d). The 
shadowing standard deviation is 6dB. Thermal noise power 
spectral density is — 174dBin/ Hz, the system bandwidth 
w is 20MHz. The minimal data rate R m in is 20Mbits / s, 
corresponding to a minimal SINR of OdB. 

We first consider an hexagonal network with 12 BSs and 
inter-site distance of 500 m, using a wrap-around to avoid 
border effects and ensure that all BSs are symmetric. Namely 
the BSs are placed on a torus. We are interested in stability of 
the operating point P* , P* = AQdBm , 1 < i < I. The corre- 
sponding target value for Gj is Gj = 80%, 1 < i < I. Without 
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coordination, this operating point is unstable as shown by 
calculating the Jacobian of G at P* using finite differences and 
calculating its eigenvalues. Another illustration of instability is 
given by plotting trajectories of the corresponding (non-linear) 
ODE starting in a neighbourhood of P*. On figures [2] and [3] 
we represent the transmitted powers and coverage probabilities 
respectively as a function of time obtained by discretization 
of the ODE. Only BSs 1, 2 and 3 are represented for clarity 
, where BSs 2 and 3 are neighbors of BS 1. Clearly, P* is 
not stable and the plotted solution does not remain close to 
P*. On figures |4] and [5] we show the same quantities when the 
coordination mechanism is applied with C = — (JG(P*)) T 
as explained in the previous section. Indeed, the solution stays 
close to P* at all times. 

Practical networks do not exactly follow an hexagonal 
model due to non-uniformity of traffic and scarcity of sites 
with good propagation characteristics. A popular model to 
take into account this irregularity is to assume that BSs 
locations follow a Poisson point process on the plane. Based 
on measurements from operational networks, results in lfl4l 
suggest that from the point of view of coverage, the Poisson 
model is pessimistic while the hexagonal model is optimist 
and reality lies somewhere in-between. We show that, like 
in the hexagonal case, instability occurs with a non-negligible 
probability. Hence instability is not an artifact of the hexagonal 
model. 

We use the following procedure. For each snapshot we 
generate BS locations according to a Poisson process on a 
square area of Akra 2 , and we find a point P* = {-P*} 1<? ;</ 
at which all BSs have the same coverage probability. To be 
consistent with the hexagonal model, we choose the neighbors 
of a BS to be the 6 closest BSs. We calculate the Jacobian 
matrix of G at P* to assess its stability. We simulate 100 
snapshots to estimate the probability of observing an unstable 
network. Figure [6] shows the probability of instability for 
different values of the number of BSs per unit of surface. There 
is a non-negligible probability of instability, and this probabil- 
ity rapidly goes to 1 when the network becomes denser. An 
intuitive explanation is that when the network gets denser, BS 
become closer to each other, so that the coupling between the 
corresponding SON mechanisms becomes stronger and causes 
instability. 

V. Conclusion 

In this paper we have studied the problem of coordinating 
multiple SON entities operating in parallel. Using tools from 
control theory and Lyapunov stability, we have proposed a co- 
ordination mechanism to stabilize the system. The mechanism 
can be implemented in a distributed fashion so it is scalable 
with respect to the number of SONs. We have shown that 
the mechanism remains valid in the presence of measurement 
noise, using stochastic approximation. Instability and coordi- 
nation in the context of wireless networks have been illustrated 
with two examples. We have shown that even for two control 
loops, instability can occur, and the influence of network 
geometry has been investigated. An interesting continuation 
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Fig. 2. Hexagonal network, no coordination, transmitted powers as a function 
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of this work would be to investigate coordination without 
linearity/linearization. Coordination for non-linear systems is 
more challenging because we cannot rely on local analysis and 
Lyapunov functions which are quadratic forms. 
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Appendix A 
Asymptotic behavior of ODEs 

Consider the ODE: 



6 = F(6), 



(34) 



which we assume to have a unique solution for each initial 
condition defined on R + . We write $(t,9(0)) the value 
at t of the solution for initial condition 9(0). We denote 
by du(9) = inf \\6 — u\\ the distance to set U. We say 

that U is invarient if 9(0) £ U implies <£>(£, 9(0)) £ U 
, t £ K + . We say that U is Lyapunov stable if for all 
8\ > there exists 82 > such that du(9(0)) < 82 implies 
du(^(t,9(0))) < <5i , t € M + . A compact invariant set U 
is an attractor if there is an open invariant set A such that 
9(0) e A implies d u (Mt, 6(0))) . A is called 

t— S- + OG 

the basin of attraction. A compact invariant set U is locally 
asymptotically stable if it is both Lyapunov stable and an 
attractor. If its basin of attraction A is equal to the whole 
space then U is globally asymptotically stable. Asymptotic 
stability is often characterized using Lyapunov functions. A 
positive, differentiable function V : M. 1 — > R, is said to be 
a Lyapunov function if t M> V($(t, 9(0))) is decreasing, and 
strictly decreasing whenever V(§(t, 9(0))) > 0. Then the set 
of zeros of V is locally asymptotically stable. If we add the 



condition V(9) — > +00, then we have global asymptotic 

||0||-H-oo 

stability. 

Appendix B 
Linear ODEs 

Consider the ODE: 

9 = A{6-9*). (35) 
Its solution has the form: 

9{t)=6* + e tA {9{{))~9*). (36) 

We denote by -< positive negativity for symmetric matrices. 
9* is asymptotically stable iif all the eigenvalues of A have 
a strictly negative real part. Alternatively, asymptotic stability 
applies iif there exists ~< X such that A T X + XA -< 0. In 
this case, V(9) = (9 — 9*) T X(9 — 9*) is a Lyapunov function 
for the ODE. The reader can refer to Q for the linear matrix 
inequality approach to stability. 



